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ABSTRACT 

We systematically analyze the prompt emission and the early afterglow data 
of a sample of 31 GRBs detected by Swift before September 2005, and estimate 
the GRB radiative efficiency BAT's narrow band inhibits a precise determination 
of the GRB spectral parameters, and we have developed a method to estimate 
these parameters with the hardness ratio information. The shallow decay com- 
ponent commonly existing in early X-ray afterglows, if interpreted as continuous 
energy injection in the external shock, suggests that the GRB efficiency previ- 
ously derived from the late-time X-ray data were not reliable. We calculate two 
radiative efficiencies using the afterglow kinetic energy Ek derived at the puta- 
tive deceleration time (idee) and at the break time (i&) when the energy injection 
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phase ends, respectively. At t b XRFs appear to be less efficient than normal 
GRBs. However, when we analyze the data at t^ ec XRFs are found to be as 
efficient as GRBs. Short GRBs have similar radiative efficiencies to long GRBs 
despite of their different progenitors. Twenty-two bursts in the sample are iden- 
tified to have the afterglow cooling frequency below the X-ray band. Assuming 
e e = 0.1, we find r/^itb) usually < 10% and rj^itdec) varying from a few percents 
to > 90%. Nine GRBs in the sample have the afterglow cooling frequency above 
the X-ray band for a very long time. This suggests a very small e# and/or a very 
low ambient density n. 

Subject headings: gamma rays: bursts - shock waves - radiation mechanisms: 
nonthermal- method: statistics 



1. Introduction 

Gamma-ray bursts (GRBs) are believed to be the most luminous electromagnetic ex- 
plosions in the universe. These erratic, transient events in 7-rays are followed by long-lived, 
decaying afterglows in longer wavelengths. The widely accepted model of this phenomenon 
is the fireball model (Meszaros 2002; Zhang & Meszaros 2004; Piran 2005), which depicts 
the observed prompt gamma-ray emission as the synchrotron emission from the internal 
shocks in an erratic, unsteady, relativistic fireball (Rees & Meszaros 1994), and interprets 
the broadband afterglow emission as the synchrotron emission from an external shock that 
expands into the circumburst medium (Meszaros & Rees 1997; Sari et al. 1998). The GRB 
radiative efficiency, which is defined as 

% - (i) 

is an essential quantity to understand the nature of the bursts. Here E 1 is the isotropic 
gamma-ray energy 1 and Ex is the isotropic kinetic energy of the fireball right after the 
prompt gamma-ray emission is over. It gives a direct measure of how efficient the burster 
dissipates the total energy into radiation during the GRB prompt emission phase. 

In the pre- Swift era, only the late time fireball kinetic energy Ek was derived or es- 
timated using the late time afterglow data. Two methods have been proposed. The most 



1 Notice that the notation E 1 is different from that used in some other papers (e.g. Frail et al. 2001) that 
denotes the geometry-corrected gamma-ray energy. 
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adequate one is through broadband afterglow modeling (e.g. Panaitescu & Kumar 2001). 
This method requires well-sampled multi-wavelength afterglow data. The method thus can 
only be applied to a small sample of GRBs. A more convenient method is to use the X-ray 
afterglow data alone (Freedman & Waxman 2001; Berger et al. 2003; Lloyd-Ronning & 
Zhang 2004). At a late enough epoch (e.g. 10 hours after the burst trigger), the X-ray band 
is likely above the cooling frequency, so that the X-ray flux gives a good measure of the de- 
generate quantity e e E K , where e e is the fraction of the electron energy in the internal energy 
of the shock. If e e could be estimated, one can then derive E K directly from the X-ray data. 
Assuming a simple extrapolation of the late time light curve to earlier epochs, Lloyd-Ronning 
& Zhang (2004) took into account the fireball radiative loss correction to estimate E K right 
after the prompt emission phase, and estimated the efficiency of 17 GRBs/XRFs observed 
in the pre- Swift era. They discovered a shallow positive correlation between i] y and E 1 or 
E p . According to this shallow correlation, softer, under-luminous bursts (e.g. X-ray flashes) 
tend to have a lower radiative efficiency Similar conclusions were also drawn by Lamb et al 
(2005). 

These previous GRB efficiency studies employing the late afterglow data inevitably in- 
troduce some uncertainties on the Ex measurements, including the possible corrections of 
radiative fireball energy loss and additional energy injection in the early afterglow phase. 
In order to reduce these uncertainties very early afterglow observations are needed. The 
successful operation of NASA's Swift GRB mission (Gehrels et al. 2004) makes this possible. 
Very early X-ray afterglow data for a large number of GRBs have been recorded by the 
X-Ray Telescope (XRT) on board Swift (Burrows et al. 2005a). A large sample of early 
X-ray afterglow data have been collected, typically 100 seconds after the triggers. These ob- 
servations indeed show novel, unexpected behaviors in the early afterglow phase (Tagliaferri 
et al. 2005; Burrows et al. 2005b; Nousek et al. 2006; Zhang et al. 2006; O'Brien et al. 
2006), and make it feasible to more robustly estimate E K and hence, ?? 7 . 

The XRT light curves of many bursts could be synthesized to a canonical light curve that 
is composed of 5 components (Zhang et al. 2006), an early rapid decay component consistent 
with the tail of the prompt emission, a frequently-seen shallow-than-normal decay component 
likely originated from a refreshed external forward shock, a normal decay component due to 
a freely expanding fireball, an occasionally-seen post jet break segment, as well as erratic 
X-ray flares harboring in nearly half of Swift GRBs that are likely due to reactivation of the 
GRB central engine. The most relevant segment for the efficiency problem is the shallow 
decay component (Zhang et al. 2006; Nousek et al. 2006). Most of the Swift/XRT afterglow 
light curves in our sample have such an early shallow decay segment. The origin of this 
segment is currently not identified. If it is due to continuous energy injection, the initial 
afterglow energy Ek must be significantly smaller than estimated using the late time data. 
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The previous efficiency analysis using late X-ray afterglow data then tend to overestimate 
E K , and hence, underestimate ?7 7 . The early tight UVOT upper limits for many Swift GRBs 
are also likely related to this shallow decay component (Roming et al. 2006). It is therefore 
of great interest to revisit the efficiency problem using the very early XRT data. 

X-ray flashes (XRFs, Heise et al. 2001; Kippen et al. 2002) naturally extend long- 
duration GRBs into the softer and fainter regime (e.g. Lamb et al. 2005; Sakamoto et al. 
2006). It is now known that the softness of the bursts is not due to their possible high 
redshifts (Soderberg et al. 2004, 2005). The remaining possibilities include either extrinsic 
(e.g. viewed at different viewing angles, Yamazaki et al. 2002, 2004; Zhang et al. 2004a,b; 
Liang & Dai 2004a; Huang et al. 2004), or intrinsic (e.g. different burst parameters such as 
Lorentz factor, luminosity, etc, Dermer et al. 1999; Kobayashi et al. 2002; Meszaros et al. 
2002; Zhang & Meszaros 2002c; Huang et al. 2002; Rees & Meszaros 2005; Barraud et al. 
2005) reasons. The radiative efficiency of XRFs may provide a clue to identify the correct 
mechanism. Previous analyses of late time X-ray data suggest that the XRFs typically have 
lower radiative efficiencies (Soderberg et al. 2004; Lloyd-Ronning & Zhang 2004; Lamb et 
al. 2005). It is desirable to investigate whether this is still true with the early afterglow 
data. 

Recently, afterglows of several short-duration GRBs have been detected (Gehrels et al. 
2005; Fox et al. 2005; Villasenor et al. 2005; Hjorth et al. 2005; Barthelmy et al. 2005a; 
Berger et al. 2005c). The data suggest that they are distinct from long GRBs and very likely 
have different progenitor systems. It is therefore of great interest to explore the radiative 
efficiencies of short-hard GRBs (SHGs) and compare them with those of long GRBs. 

In this paper we systematically analyze the prompt emission and the early afterglow 
data for a sample of 31 GRBs detected by Swift before September 2005 through reducing the 
Burst Alert Telescope (BAT, Barthelmy et al. 2005b) and XRT data. We present the sample 
and the BAT/XRT data analysis methods in §2. The gamma-to-X fluence ratio (R^/x) is 
an observation-defined apparent GRB efficiency indicator. We perform statistical analysis of 
this parameter in §3. In §4, we perform more detailed theoretical modeling to estimate E K 
(§4.1) and r^ 7 (§4.2). Our results are summarized in §5 with some discussion. Throughout 
the paper the cosmological parameters H = 71 km s _1 Mpc -1 , VLm = 0.3, and = 0.7 
have been adopted. 
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2. Data 

Our sample includes 31 GRBs observed by Swift before September 1, 2005. The BAT- 
XRT joint early X-ray afterglow light curves of these bursts and the detailed data reduction 
procedures have been presented in O'Brien et al. (2006). 

2.1. Prompt Gamma-rays 

The BAT data have been processed using the standard BAT analysis software (Swift 
software v. 2.0). The BAT-band gamma-ray fluence S lj0 b s could be directly derived from 
the data. For the purpose of estimating GRB efficiency, on the other hand, one needs to 
estimate the total energy output of the GRB, which requires to extrapolate S 7j0 b s to a wider 
bandpass to get S* 7 (1 — 10 4 keV adopted in this paper). This requires the knowledge of the 
spectral parameters of the prompt emission. 

It is well known that the GRB spectrum is typically fitted by a Band function (Band et 
al. 1993), which is a smoothly-joint-broken power law characterized by two photon indices 
Ti and T 2 (with the convention N(E)dE oc E r dE adopted throughout the text) and a break 
energy E . The peak energy of the vf u spectrum is E p = (2 + Ti)E . In order to derive these 
parameters, the observed spectrum of a burst should cover the energy band around E Q . The 
spectra of both long and short GRBs observed by CGRO/BATSE (covering 20-2000 keV) 
are well fitted by the Band-function, with the typical values of Ti ~ —1, T 2 ~ —2.3, and 
E p ~ 250 keV (Preece et al. 2000; Ghirlanda et al. 2003). The spectra of XRFs could be 
also fitted by the Band-function, typically with l?i ~ — 1 and a lower T 2 than typical GRBs 
(Lamb et al. 2005, Sakamoto et al. 2005, 2006; Cui et al. 2005). BAT has a narrower 
energy band (15-150 keV) than BATSE and HETE-2. The typical E p of the bright BATSE 
sample is well above the BAT band. BAT's observations therefore cannot well constrain E p 
and T 2 for many GRBs. Most observed spectra by BAT are well fitted by a simple power 
law. However, this is due to the intrinsic limitation of the instrument. Four GRBs in our 
sample, i.e. GRBs 050401, 050525A, 050713A, and 050717, were simultaneously observed by 
Konus-Wind (with an energy band of 20-2000 keV). The spectra of these bursts could be also 
well fitted by a Band-function or a cutoff power law spectrum (Golenetskii et al. 2005a-d, 
Krimm et al. 2006). We assume that the broad band spectra of all the bursts in our sample 
could be fitted by the Band-function and then make the corrections to the observed fluences. 
Our procedure to derive spectral parameters are as follows. For more detailed data analysis 
procedure and burst parameters, we refer the readers to the Swift first catalog paper (L. 
Angelini et al. 2006, in preparation). 
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We first fit an observed spectrum with a Band-function, a power law with exponential 
cutoff, and a simple power law, respectively. By comparing reduced x 2 °f these fits, we 
pick up the best-fit model among the three. In our sample, only 4 GRBs (050128, 050219A, 
050525A and 050716) could be well fitted by a Band function, if one assigns T 2 in the range of 
—5 to — 2 2 . Due to the great uncertainty of T 2 , the cutoff power law model could also fit these 
4 bursts, with the cutoff energy being in the BAT band. The rest 27 bursts in the sample are 
best fitted by a simple power law, with the photon index Y PL ranging from ~ — 1 to ~ — 3 
(see Table 1). For the 8 GRBs with the Band-function parameters available (GRBs 050128, 
050219A, 050401, 050525A, 050713A, 050716, 050717, and 051221) either from a fit to the 
BAT data or from a joint BAT- Konus- Wind fit, we make straightforward extrapolation of 
Sj,obs to derive S* 7 in the 1 — 10 4 keV band. For the rest 24 bursts whose observed spectra are 
fitted by a simple power law, the extrapolation is not straightforward. Generally we employ 
the hardness ratio information to place additional constraints to the spectral parameters. 
The analyses are carried out on case-to-case base. Nonetheless, depending on the value of 
T PL , one could crudely group the cases into three categories. 

Case I: —2.3 < T PL < —1.2 (16 out of 24 in the sample). Since for typical bursts in the 
BATSE sample one has Ti ~ —1 and T 2 ~ —2.3 (Preece et al. 2000), it is expected that 
a rough fit to the Band function by a simple power law would lead to Y PL in this range. 
The break energy E of these bursts should be within or near the edges of the BAT band. 
The hardness ratio (HR) , which in our analysis is defined as the ratio of the fluence in the 
50-100 keV band to that in the 25-50 keV band, could be directly measured from the simple 
PL fit model. Theoretically, on the other hand, HR is a function of Ti, T 2 and E for the 
Band-function model (Cui et al. 2005, see Fig.l). One can then in principle apply HR as 
another agent to constrain the spectral parameters by requiring 

HR mod = HR obs , (2) 

where H R mod and HR obs are the hardness ratio derived from the Band-function model and 
from the data, respectively. To proceed, we first assign a set of "standard" guess values to 
the spectral parameters, e.g. Ti = —1, T 2 = —2.3, and E p = 250 keV, and then perform 
a Band-function fit to the data. We then derive T\, T 2 , and Eq from the best fit, which 
usually deviate from the guess values. Generally, these parameters have very large error 
bars. We then apply the HR criterion to the results. If the calculated HR mod using the 
best-fit parameters match HR obs well, this set of parameters are taken. Otherwise, we 
adjust spectral parameters to achieve the best match. The process is eased thanks to several 



2 Although T 2 2.3 for typical GRBs, T 2 could be as low as 5 for XRFs (e.g. Cui ct al. 2005; 

Sakamoto et al. 2005). 
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Fig. 1. — The HR — E relation for a Band-spectrum. Different sets of (ri, T 2 ) have been 
plotted. 
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properties of the HR — E relation (Fig.l). First, when E p is high enough (say, higher than 
30 keV), HR is essentially independent on T 2 . Also for T 2 < —2, most energy is emitted 
around E p and the extrapolated broad band fluence is insensitive to T 2 . We therefore take 
the best-fit T 2 value or fix it to ~ —2.3 when T 2 is poorly constrained. HR mod therefore 
mainly depends on T\ and Eq. Another interesting feature is that as one decreases I\ (i.e. 
softer spectrum), the corresponding E would increase given the same observed spectrum. 
The resulting E p = (2 + Ti)E , on the other hand, is not very sensitive to r x since the 
variations of r\ and E cancel out each other. As a result, we simply adjust r x to re-fit 
the spectrum until HR mod is consistent with HR obs within the error range. These spectral 
parameters (r 1; E p , T 2 ) are then taken to perform extrapolation to estimate S 7 . Since the 
determinations of both Ti and T 2 are not fully based on fitting procedures, it is very difficult 
to quantify their errors, and we only report their estimated values. The error of E is taken 
whenever possible based on the best Band-function fit (see Table 1). 

Case II: T PL > -1.2 (2 out of 24 in the sample: GRBs 050726 and 050826). In this 
case, Eq should be far beyond the BAT band, and BAT only covers the low energy part of 
the spectrum (E < E ). One has Ti ~ T PL or slightly larger (if E p is not very far above the 
band). It is very difficult to estimate E (E p ) in this case. Nonetheless, one could use HR obs 
to pose some constraints. Figure 1 suggests that HR mod converges to a maximum value at 
high E p s given a certain T 1 . We first let T 1 = T PL and check whether HR obs is consistent 
with HR mod . For both bursts in our sample one has HR obs > HR™° d . This suggests that 
Ti should be larger than T PL . We then gradually increase Ti until HR™° d (Ti) becomes 
consistent with HR obs . Using this T 1 one can then fit for E (and E p ). According to Fig.l, 
there is great degeneracy of E (E p ) at HR™° d . The fitted E (E p ) therefore has very large 
errors and is unstable. In practice, one could only set a lower limit of Eq (E p ) below which 
HR mod starts to deviate from HR obs . As a result, only lower limits of E p of these two bursts 
are reported in Table 1. 

Case III: T PL < —2.3 (6 out of 24 in the sample). In this case, the burst is likely an XRF 
with E p near or below the low energy end of BAT. BAT's observation likely only covers the 
high energy part of the spectrum (E > E ). In order to constrain the spectral parameters 
using the hardness ratio data, we assume T 1 = — 1 and T 2 = Y PL and fit for E by requiring 
HR mod (E ) = HR obs . According to Fig.l, near HB™&, HR mod is rather insensitive to E 
(E p ). For most of the cases, one could find a solution of E , but with very large errors. The 
solutions are also unstable. For these cases, we do not report errors in Table 1, but only 
report the best fit value with a ~ symbol 3 . For two cases (GRBs 050416A and 050819), 



3 Similarly, GRBs 050319, 050724, 050801 in the Case I also have unstable solutions, so that the errors of 
their E p 's are not reported. 
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there is no solution since HR obs < HR™°%. Similar to the Case II, we then lower T 2 until a 
solution is found. Due to the degeneracy of E (E p ) with HR, only upper limits of E Q (E p ) 
are found, which are reported in Table 1. 

With the spectral parameters derived from the above method, we have extrapolated the 
observed 15-150 keV fluence S Jj0 b s to a broader energy range (1 — 10 4 keV). The derived S y 
is regarded as the total energy output during the prompt phase for further efficiency studies. 
The results of prompt emission data are reported in Table 1. We caution that due to the 
intrinsic instrumental limitation, the uncertainties of the results are large. Nonetheless, we 
have made the best use of the available data (especially the hardness ratio information) to 
derive the parameters. Figure 2 displays the robustness of the method. Fig. 2a shows the 
criterion adopted in analyzing each burst (eq. [2]). Fig.2b shows how the E p s derived with 
our method compares with the E£ %t derived from the Konus-Wind or HETE-2 data for 8 
bursts. The result suggests that our derived E p meets E^ lt well for moderate E p s (i.e. those 
falling into the BAT band), but deviate from E^ lt when E p is very large. Notice that the 
reported E p errors in our method are derived from the best fit by fixing T\ and T2- The real 
errors should also include the uncertainties of Ti and T2- In such cases, the errors of E p s 
could be larger to be more consistent with E^ lt in the high-i? p regime. However, due to the 
difficulty of estimating these errors, they are not included in reported errors. 

In Fig. 3, we show the distributions of HR and E p of our sample. For the E p distribution 
(right panel), those of the BATSE sample (Preece et al. 2000) and the HETE-2 sample (Lamb 
et al. 2005) are also plotted for comparison. Our sample is generally consistent with the 
HETE-2 sample, and tends to be softer than the BATSE sample (as is expected because of 
a softer bandpass of BAT as compared with BATSE). If one defines XRFs as those bursts 
with HR obs < 1, the number ratio of the XRFs and the GRBs in our sample is ~ 1 : 2, being 
consistent the HETE-2 result (Lamb et al. 2005). An interesting feature is the marginal 
bimodal distribution of XRFs and GRBs, which is consistent with the previous result derived 
with the HETE-2 data (Liang k Dai 2004, cf. Sakamoto et al. 2005). In Fig.4, we plot the 
distribution of as compared with the BATSE results. The energy range of S 1 is re-scaled 
to 20-2000 keV (BATSE's energy band). We find that the two distributions are generally 
consistent with each other, except that Swift GRB sample extends the BATSE sample to 
lower fluences. This is expected because Swift is more sensitive than BATSE. 
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Fig. 2. — (a) Comparison of HR mod with HR obs . (b) Comparison of E p derived from our 
method with E^ lt derived from the joint fit using BAT, Konus-Wind, or HETE-2 data for 
eight bursts, including GRBs 050128, 050219A, 050401, 050525A, 050713A, 050716, 050717, 
and 051221. The solid line is E p = 



-i 1 1 1 1 1 r 



i — 1 — r 



1.0 - 



0.8 - 



0.6 - 



0.4 - 



0.2 - 



0.0 L-L- 1 ' 1 ' 1 ' 1 ' L- 1 - 

-0.4 -0.2 0.0 0.2 0.4 



log Hft 



bs 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



log EJ keV 



Fig. 3. — The distributions of HR obs (left) and E p (right) derived from our method. For 
the E p distribution, the histograms of BATSE (dotted line) and HETE-2 (step-dashed line) 
samples are also plotted. 
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Fig. 4. — Comparison of log 5^ distribution (in the 20-2000 keV band) of our sample (solid 
histogram) to that of the BATSE sample (dashed histogram). 



-13- 



2.2. X-ray afterglows 

The X-ray afterglow light curves and the spectra of the bursts in our sample have 
been presented by O'Brien et al. (2006). Among the 5 components of the synthetic X-ray 
lightcurve (e.g. Zhang et al. 2006), the steep decay component is due to the GRB tail 
emission (Kumar & Panaitescu 2000; Zhang et al. 2006) and the X-ray flares are due to the 
late central engine activity (Burrows et al. 2005b; Zhang et al. 2006; Fan & Wei 2005; Ioka 
et al. 2005; Liang et al. 2006; Wu et al. 2006; King et al. 2005; Perna et al. 2006; Proga & 
Zhang 2006; Dai et al. 2006). We therefore identify the steep decay component as well as 
the X-ray flares and remove them from the light curve contribution. We then fit the light 
curves by either a broken power law for those bursts with shallow-to-normal transition, or 
by a single power law otherwise. Table 2 lists the X-ray data and the fitting results of our 
sample. 

3. Prompt Gamma-ray fluence vs. X-ray Afterglow Fluence 

In order to calculate the absolute value of the GRB radiative efficiency ^ 7 (eq.[l]), both 
Ej and E^ need to be derived. This requires detailed modeling of Ex and the redshift infor- 
mation. The results depend on some unknown parameters (e.g. the shock electron/magnetic 
field equipartition factors, e e , e#, etc), which we discuss in detail in §4. Nonetheless, using 
the directly measured quantities listed in Tables 1 & 2, one can analyze the relative ener- 
getics between the prompt emission and the afterglow. The prompt emission fluence S y is 
a rough measure of the prompt emission energetics. According to the standard afterglow 
model and assuming that electron spectral index is p ~ 2, that the X-ray band is above both 
the typical synchrotron emission frequency v m and the synchrotron cooling frequency v c , 
and that the inverse Compton cooling is unimportant, the afterglow kinetic energy could be 
roughly indicated by the quantity e~ 1 F l/ x(t)t, where t is a particular epoch in the afterglow 
phase (e.g. Freedman & Waxman 2001; Berger et al. 2003; Lloyd- Ronning & Zhang 2004). 
Since the quantity Sx(t) = F Vt x(t)t also has the dimension of fluence, the S^-to-Sxif) ratio 
could give a rough indication of the relative energetics between the prompt emission and the 
afterglow. The unknown redshift essentially does not enter the problem. 

The shallow decay phase commonly observed in Swift GRBs has been generally inter- 
preted as a refreshed external shock (e.g. Rees & Meszaros 1998; Dai & Lu 1998; Panaitescu 
et al. 1998; Kumar k Piran 2000; Sari & Meszaros 2000; Zhang & Meszaros 2001, 2002a; 
Dai 2004; Zhang et al. 2006; Nousek et al. 2006; Panaitescu et al. 2006; Granot & Kumar 
2006). Within this interpretation, the kinetic energy of the afterglow E K increases with time 
for an extended period. This brings extra complication to the efficiency problem. One needs 
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to identify at which epoch the corresponding Ek represents the kinetic energy left over right 
after the prompt gamma-ray emission. This is a model-dependent problem, and we take an 
approach to accommodate different possibilities. We pay special attention to two epochs. 
One is the break time t& at the shallow-to-normal-decay transition epoch, which corresponds 
to the epoch when the putative injection phase is cover. Within the injection interpretation, 
another important time is the fireball deceleration time (tdcc), which is usually earlier or 
around the first data point in the shallow decay phase. The kinetic energy at tdcc [Exit dec)] 
is relevant to the GRB efficiency problem, if the injected energy during the shallow decay 
phase is due to a long-term central engine (e.g. Dai k Lu 1998; Zhang k Meszaros 2001; Dai 
2004), since the bulk of kinetic energy is injected at later epochs. In the scenario that the 
injection is due to an instantaneous injection with variable Lorentz factors (Rees k Meszaros 
1998; Kumar k Piran 2000; Sari k Meszaros 2000; Zhang k Meszaros 2002a; Granot k Ku- 
mar 2006), the total kinetic energy of the outflow right after the prompt emission is over 
should be defined by the kinetic energy measured at tj, [Exitb)] when the injection phase is 
over. However, the kinetic energy of the ejecta that gives rise to the gamma-ray emission may 
be still roughly E K (td ec ), since the gamma-ray emission from the low-Lorentz-factor ejecta 
can not escape due to the well-known compactness problem (e.g. Piran 1999). Nonetheless, 
some other scenarios do not interpret the shallow decay as additional energy injection (e.g. 
Eichler k Granot 2006; Toma et al. 2006; Ioka et al. 2006; Kobayashi k Zhang 2006). In 
some of these cases, Exitb) is the more relevant quantity to define the GRB efficiency, e.g. 
in the models invoking precursors (e.g. Ioka et al. 2006) and the models involving delayed 
energy transfer (e.g. Kobayashi k Zhang 2006; Zhang k Kobayashi 2005). In this paper we 
use both Ex{tdec) and Ex{tb) to define GRB radiative efficiencies. 

The injection break time t& can be directly measured from the light curves. The fireball 
deceleration time (for burst durations shorter than this time scale, the so-called thin shell 
regime) tdcc ~ 5(1 + ,2)(-Efr- i 52 / M) 1 / 3 (7o/300) -8 / 3 , on the other hand, is not directly measured 
and is very likely buried beneath the steep-decay prompt emission tail component. Here 
Ek$2 = E K /10 52 ergs (the convention Q n = Q/10 n in cgs units is adopted throughout the 
paper), n is the ambient medium density in unit of 1 proton / cm 3 , and 70 is the initial Lorentz 
factor of the fireball. Without knowing 7 (noticing the sensitive dependence on this unknown 
parameter), one cannot accurately estimate this time with the observables. For typical 
parameters (e.g. z — 1, 70 = 150, and Ex^/n = 1) we get td ec ~ 60s. Considering also the 
thick shell regime (i.e. the burst duration is longer than the above critical time, Kobayashi 
et al. 1999), we finally roughly estimate the deceleration time as td cc ~ max(60 s, T 90 ). 

In Figure 5 we plot <Sx(tdec) ( re d dots) and Sx(tb) (black triangles) against S T The large 
differences between Sx(tdec) an d Sx{h) indicate that significant energy injection happens in 
many bursts. One interesting signature evidenced in Fig. 5 is that S 7 is positively correlated 
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to Sx- This is consistent with the previous knowledge that the radiated energy is positively 
correlated to the afterglow kinetic energy. For the case of Sx{tb), the S 7 — Sx correlation 
slope is 0.78 ± 0.17, which is shallower than unity. This suggests that considering the total 
kinetic energy when the energy injection phase is over [Ex{tb)\, the fainter bursts are not 
as efficient as brighter ones in converting kinetic energy into radiation. This is consistent 
with the pre-Swift finding of a shallow correlation between the efficiency r? 7 and E 1 (Lloyd- 
Ronning & Zhang 2004; Lamb et al. 2005). Since generally there is a positive correlation 
between the and E p (Lloyd et al. 2000; Amati et al. 2002; Liang et al. 2004a; Lamb et 
al. 2005), this also suggests that softer bursts (e.g. XRFs) tend to have lower efficiencies 
(Soderberg et al. 2004; Lloyd- Ronning & Zhang 2004). For the case of t dec , however, S 1 — S x 
correlation slope is very close to unity (1.17 ± 0.22). This means that the efficiency defined 
by Exitdec) (the "true" efficiency in the models invoking additional energy injection) is not 
sensitively related to E y , and hence, to E p according to the positive E 1 — E p correlations. 
This means that XRFs are not intrinsically inefficient GRBs if the injection hypothesis is 
true. 

To show the effect more clearly, we plot the ratios R(tdec) = S 1 / Sx{tdec) (red dots) 
and R{tt b ) = S^/Sxitb) (black triangles) against the hardness ratio HR obs (Fig.6). It is 
found that while the soft bursts (e.g. XRFs) tend to have a lower gamma-to-X ratio (and 
hence, lower radiative efficiency) than the typical GRBs at t b (a shallow linear dependence 
with index 0.61 ± 0.21, they do not differ too much from hard GRBs at tdcc- This means 
that XRFs are radiatively as efficient as hard GRBs. Such a possibility has been speculated 
by Lloyd-Ronning & Zhang (2004) and was first recognized by Schady et al. (2006) when 
analyzing the early UVOT data of XRF 050406. Now we extend the analysis to a larger 
sample and verify that it is common for other XRFs as well. Such a conclusion is strengthened 
by a more careful treatment of radiative efficiency in §4.2, and we discuss the implications 
of this result for the XRF models in §5. 
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Fig. 5. — The extrapolated prompt emission fluence S~, against the X-ray fluence at td ec 
(dots) and at t b (triangles). The solid and dashed lines, with slopes of 0.78 and 1.17, are 
the best fit to the data at t b and ^ ec , respectively. The short-hard GRBs (SHGs) and X-ray 
flashes (XRFs) are marked. 
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Fig. 6. — Ratios of S^/Sx against the hardness ratio HR. The circles represent the data at 
tdeci and the triangles represent the data at £&. The short-hard GRBs (larger solid circles) 
and the XRFs (open circles) at td ec are marked. The straight line is the best fit of the data 
at tb, which shows a shallow correction between the two quantities, with a linear correlation 
coefficient 0.61 and a chance probability < 0.01. 
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4. GRB efficiency 

4.1. Theoretical Models 

In this section we explicitly derive the radiative efficiency for the GRBs in our sample ac- 
cording to eq.(l). The isotropic prompt emission energy E 1 is derived from the extrapolated 
1-10000 keV band ffuence (S' 7 ) according to 

E 1 = 47^^(1 + ^-1 

= 1.3 xlO 51 erg D 2 L28 (1 + z)~ 1 S 1 ,_ 6 (3) 

where Dl is the luminosity distance. 

The derivation of Ek requires detailed afterglow modeling. Regardless of whether there 
is indeed a long-lasting central engine, the energy injection process could be mimicked by 
introducing an "effective" long lasting central engine with luminosity L = L (t/t )~ q . For 
the varying Lorentz factor injection model the Lorentz factor index s could be related to 
an effective q through q = (10 — 2s)/ (7 + s) for an ISM medium and q = 4/(3 + s) for a 
wind medium (Zhang et al. 2006). The injection then results in an evolving kinetic energy 
E K oc t {l ~ q \ After the energy injection is over, the fireball can be described by the standard 
afterglow model. At any epoch during the injection phase, the afterglow emission level could 
be calculated by taking the kinetic energy at that time. So in our treatment, we still use the 
standard afterglow model to derive various parameters as functions of E K , bearing in mind 
that Ek may be time-dependent. 

For a constant density medium, the typical synchrotron emission frequency, the cooling 
frequency and the peak spectral flux read (Sari et al. 1998; coefficients taken from Yost et 
al. 2003) 4 

u m = 3.3xl0 12 Hz(^y ( l + ,)^ (4) 

v c = 6.3 x 10 15 Hz(l + z)- l ' 2 (l + Y)- 2 ef / J 2 E K ]lln- 1 r d 1/2 (5) 
F„ imax = 1.6 mJy(l + z)D 2 2 e 1 J 2 _ 2 E K>52 n 1 ^ . (6) 

4 Lloyd-Ronning & Zhang (2004) adopted the coefficients from Hurley et al. (2002). The v m coefficient is 
larger than adopted here. This effect, together with the ignorance of the Inverse Compton (IC) effect, leads 
to systematic underestimating of Ek and overestimating of ry 7 , as also pointed out by Fan & Piran (2006) 
and Granot et al. (2006). This systematic deviation does not affect the global dependences of 77 7 on other 
parameters, as are reproduced in this paper. 
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Here td is the observer's time in unit of days, 

y=[-l + (l+4 We AB) 1/2 ]/2 (7) 

is the IC parameter, where r/i = min[l, (z/ c / ' v m ) <y2 ~ p ^ 2 ] (Sari & Esin 2001), and r/ 2 < 1 is a 
correction factor introduced by the Klein-Nishina correction. The latter effect was treated 
in detail by Fan & Piran (2006). Here we adopt an alternative, approximate treatment. For 
the electron Lorentz factor 7 e x corresponding to the X-ray band emission, the synchrotron 
self-IC effect is significantly suppressed in the Klein-Nishina regime for the photons with 
energy v > z/kn, where 

z/kn = rTm e c 2 7^(l + z)' 1 

- 2.4 x 10 15 Hz (1 + zY'M^t-^ 1 ^' 2 , (8) 

and h is the Planck's constant. Based on the vF v spectrum of the standard synchrotron 
emission model (Sari et al. 1998), one can roughly estimate r/2 = min[l, (^kn/^c)' 3 ^''' 2 ] for 
slow cooling {y m < u c ) and r/2 = min[l, (z/kn / Vm) L ^\ for fast cooling (z/ c < z/ m ), where the 
factors (^kn/^)' 3 _p " 2 and (^KN/^m) 1 ^ 2 denote the fractions of the photon energy density that 
contributes to self-IC in the X-ray band in the slow and fast cooling regimes, respectively. 

In previous analyses (e.g. Freedman & Waxman 2001; Berger et al. 2003; Lloyd- 
Ronning & Zhang 2004), the IC cooling was usually not taken into account. The inclusion 
of IC cooling modifies the X-ray afterglow light curves considerably (e.g. Wu et al. 2005), 
which also influences the derived GRB efficiency r^ 7 . Here we generally include the IC factor 
[power of (1 + Y)} in the treatment (see also Fan & Piran 2006). The results are reduced to 
the previous pure synchrotron-dominated case when 7 < 1. 

For about 2/3 cases in our sample (22 out of 31), the X-ray band temporal decay index 
and the spectral index are consistent with the spectral regime v > max(z/ m , u c ). This is the 
regime where E K is independent of n and only weakly depends on e B and p, and therefore 
an ideal regime to measure E K . One can derive the X-ray band energy flux as 5 

^=10 18 Hz) = F v ^yj 2 vt- 1)t2 v%- v)l2 

= 5.2 x 10~ 14 ergs s" 1 cnr 2 /^! + zf p+2)/A 

x (i + vrui-f 4 <-\^j 2 2)/4 tr 3p)/ v 18 (2 - )/2 , o) 

where ^ 

f p = 6.73 (j^j P (3-3 x 1Q-6) (f- 2 - 3 )/ 2 (10) 



5 Although our treatment is for a constant-density medium, eq.(9) is valid for more general cases (e.g. 
wind medium) since in this regime the flux does not depend on the medium density. 
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is a function of p, which is calculated in Figure 7. It peaks at ~ 1.72 when p ~ 2.12, and 
declines monotonically at large p values. For example, at p ~ 3, f p is only ~ 0.02. This gives 
a nearly 2 orders of magnitude variation for p = (2.01 — 3), and thus demands more careful 
treatments of individual bursts presumably having quite different p values. 

With eq.(9), one can derive E K at any time t d as 



E 



K,52 



vF v {y = 10 i8 Hz) 
5.2 x 10~ 14 ergs s _1 cm" 



4/(p+2) 



X (1 + Y) 4/(p+2) /- 4 /(P+2) e g-P)/(P+2) e 4(l-p)/(p+2)^ g 2(p-2)/(p+2) 



This could be reduced to E K ~ e e 1 S X D 2 L / (1 + z) for p ~ 2. Comparing with eq.(3), we can 
see that as far as the efficiency problem is concerned the redshift-dependence is very weak. 
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Fig. 7.— The function f p . 
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For nearly 1/3 cases in our sample (9 out of 31), the X-ray data are not consistent with 
being in the regime vx > max(i/ m , v c ) (which requires a = (3(3 — l)/2 for the F v oc t~ a v~ 13 
convention). The temporal decay slope in the normal decay phase is close to -1. This rules out 
the fast cooling case v c < vx < v m for both ISM and wind models. For slow cooling models 
(y m < vx < v c ), the temporal decay index derived from the wind model [a = (3(3 + l)/2 for 
the F v oc t~ a v~P convention, Chevalier & Li (2000)] is too steep compared with the data in 
our sample. One is then left with the only possibility, u m < u x < v c in the ISM model. In 
fact the observed a and (3 are consistent with being in this regime [a = (3/2)/?]. The derived 
X-ray band energy flux is then 



vF v (v = 10 iS Hz) 



x 



F (u /vv) ip - 1)/2 

± v,ma,x\' y m,/ "A } 
-13 



6.5 x lO - " ergs s _i cm _2 D^ 8 2 (l + z 



,(P+3)/4 



(p+l)/4 c p-l p? (p+3)/4^i/ 2+ (3-3p)/4 iy _(3-p)/2 



f e 1 

JP t B,~2 



f E 



K,52 



Vis' 



(12) 



This gives 



E 



K,52 — 



X 



vFJu 



10 18 Hz) 



6.5 x 10 13 ergs s 1 cm 2 



4/(p+3) 



-1.3(p-l)/(p+3) 



4/(p+3) <r -(p+l)/(p+3) <r 4(l-p)/(p+3) ri -2/(p+3) J/ _2(p-3)/(p+3) 



'B,-2 



-e,-l 



'^18 



(13) 



Inspecting eq.(5), one can draw the conclusion that in order to have u c > u x in the normal 
decay regime (t d ~ 1), e B must be very small (e.g. < (10~ 3 — 10 -4 )). Keeping a more or 
less constant e e ~ 0.1, the Y parameter (eq.[7]) does not increase significantly for a smaller 
€b, since the r] 2 parameter becomes much smaller due to the Klein- Nishina suppression. In 
order not to derive a unreasonably large Ek value (limited by the total energy budget of the 
progenitor system), the data also require a small ambient density (e.g. n < 0.1). 



4.2. Calculation Results 

Eqs. (11) and (13) suggest that the absolute value of E K depends on several unknown 
shock parameters. In order to calculate ?7 7 , the absolute value of E K is needed. This 
requires a detailed multi-wavelength study (e.g. Panaitescu & Kumar 2002; Yost et al. 
2003) to constrain unknown shock parameters as well. Limited by the X-ray data alone, one 
inevitably needs to make some assumptions on the unknown shock parameters. 

Our first step is to use the X-ray data (temporal index ax and spectral index (3x) to 
determine the spectral regime the burst belongs to. In all the cases, we choose the "normal" 
decay phase for the temporal index (a 2 for the broken power-law fit, or a\ for the single 
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power-law fit if a± is close to unity). This is because this segment has no contamination 
of energy injection (with unknown q parameter). Using ax and /3 X we check the spectral 
regime of the X-ray band by comparing the a — (3 relation of various models (e.g. Table 1 
of Zhang & Meszaros 2004). We find that within the error uncertainties of the data, most 
bursts could be grouped into two spectral regimes: (I) v x > max(z/ m , v c ) (20 out of 31); and 
(II) u rn < vx < v c in the ISM model (9 out of 31). Two bursts have unexpectedly large 
spectral indices, i.e. GRBs 050319 ((3 X = 2.02 ±0.47) and 050714B (0 X = 4.50 ±0.70). The 
spectrum is likely dominated by the contribution of the GRB tail emission and we assume 
that in the afterglow phase they are in the regime of v x > max(z/ m , u c ), the default case. 
GRBs 050215B (0 X ~ 0.5 ± 0.5) and 050716 (0 X = 0.33 ± 0.03) have very small spectral 
indices, and we assume that it is in the regime v m < vx < v c . After determining the spectral 
regimes, we derive p from (3 X (Px — —p/2 for u x > niax(z/ m , u c ) and (3 X = —{p — l)/2 for 
v m < v x < v cj- Given the large error bars usually associated with the spectral indices, 
whenever p < 2 and p > 3 we take p = 2.01 and p = 3, respectively. 

4.2.1. v x > max(z/ m , v c ) 

Since there are too many unknown parameters for the regime II bursts (eq.[13]), we 
first ignore them and focus on the regime I bursts, whose Ex essentially only depend on 
e e . Previous broadband fitting suggests that e e is typically around 0.1 (Wijers & Galama 
1999; Panaitescu & Kumar 2002; Yost et al. 2003; Liang et al. 2004; Wu et al. 2004). The 
value of €3 has a large scatter for previous bursts but nonetheless has a typical value of 0.01 
(e.g. Panaitescu & Kumar 2002). The existence of regime II bursts suggest that at least 
some bursts have very small e#. We therefore also consider the cases with a smaller ee, say, 
~ 10~ 4 . In any case E K is insensitive to in regime I. 

Our calculation procedure for regime I GRBs is as follows. 

1. Use the extrapolated S 1 in the 1-10000 keV band to derive E 1 according to eq.(3). 
Since r^ 7 is insensitive to z, we take a moderate redshift, z — 2, for those bursts whose 
redshifts are not directly measured. 

2. Use Px and the 0.3-10 keV band flux to calculate the monochromatic flux at 10 18 
Hz, F v (v = 10 18 Hz). 

3. Calculate Ek,52 with eq.(ll) at two epochs, tdec and tb- For each epoch we calculate 
two values: one value {E^\ 2 ) for (e e , e^) = (0.1,0.01), and another (£^52) f° r ( e e? e -e) — 
(0.1, 10~ 4 ). The Y parameter is searched self-consistently according to the method described 
in 84.1. 
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4. Use eq.(l) to derive rffl and at tdcc and U- 

Table 3 displays our calculation results for regime I GRBs. Equation (11) indicates that 
the apparent dependence on e# is weak. This is strengthened by the Klein-Nishina effect 
since Y does not increase significantly as is lowered. Comparing E K and E K (or 77* and 
77^) at a same epoch, we can see that the difference introduced by changing e B = 0.01 to 10~ 4 
is not significant. In Fig. 8 we show the Ek contour in the (e e ,e#) plane for GRB 050219A. 
It again shows the result is insensitive to €3, so that e e is the most sensitive parameter 
for determining Ek and hence, for determining r] T The existence of the regime II bursts 
(which requires low values of e#) suggests that if shock parameters are not too different 
from burst to burst, the e# value for the regime I bursts may be also low (e.g. the second 
parameter set, (e e ,e_e) = (0.1, 10~ 4 ), for the regime I calculations). This suggests slightly 
lower radiative efficiencies than previously estimated (typically taken as e# ~ 0.01), since a 
lower eB nonetheless slightly increases E K despite a very shallow dependence. 
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Fig. 8.— Contours of logE K in the (e e ,e B ) plane for GRB 050219A. 
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Inspecting the calculated r? 7 for the bursts in the spectral regime I, we find that at 
tdec, 10 out of 22 bursts have a radiative efficiency rj^ higher than 60%, sometimes even as 
high as 98% (GRB 050819). The rest have much lower efficiencies, sometimes only a few 
percent. At t^, on the other hand, the values of i] y are typically several percent or even 
lower. Those bursts with low efficiencies from the very beginning correspond to the cases 
without significant energy injection in the early phase. For illustration, in Fig.9 we present 
the combined BAT-XRT light curves in the XRT band for some GRBs having extremely 
high or low efficiencies at t dcc (see also O'Brien et al. 2006). It is evident that the XRT 
light curves of the high-?7 7 (trf ec ) GRBs have a very flat energy injection component and a 
prominent steeply decaying prompt emission tail (e.g. GRB 050315 and GRB 050714B). 
Those with low-77 7 (^), on the other hand, typically have a smooth transition from prompt 
emission to afterglow without a significant steep decay component and/or a shallow decay 
component due to energy injection (e.g. GRB 050401 and GRB 050712). 




Fig. 9. — Comparisons of the joint-BAT-XRT light curves in the XRT band for some GRBs 
having extremely high or low efficiencies at td ec - 
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In Fig. 10 we plot r]\ as a function of the hardness ratio HR. We again find that while 
there exists a shallow correlation between rfy and HR (with index 1.12 ± 0.90), such a 
correlation essentially disappears when r/ 7 at tdec is considered. This is consistent with the 
analysis in §2.2. The results suggest that if Ex{tdec) is the relevant afterglow energy left over 
after the prompt emission, XRFs and GRBs intrinsically have similar radiative efficiencies, 
as opposed to the conclusion drawn using late-time X-ray data only. 
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Fig. 10. — The radiative efficiency f]\ against the hardness ratio HR. The circles represent 
the data at td ec , and the triangles represent the data at i&. The short-hard GRBs (larger 
solid circles) and the XRFs (open circles) at td ec are also marked. The dashed line is the best 
fit of the data at tf,, which shows a shallow correction between the two quantities similar to 
that shown in Figure 6. 
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Two short-hard GRBs 050724 and 050813 are in the spectral regime I. The results 
in Table 3 suggest that there is no noticeable difference between short, hard GRBs and 
long, soft GRBs as far as the radiative efficiencies are concerned. The same conclusion was 
drawn for the first short GRB with X-ray afterglow detection (GRB 050509b, Gehrels et 
al. 2005; Bloom et al. 2006), and is being verified by a larger sample of short GRBs in 
accumulation. To increase the short GRB sample, we have also performed the same analysis 
to GRB 051221A (e.g. A. Parsons et al. 2006, in preparation). The results are included in 
all the Tables and figures, making the total number of short bursts in our sample to 32. The 
inclusion of this burst strengthens the conclusion that short GRBs are no different from long 
GRBs in radiative efficiencies. 



4.2.2. v m <v x < v c 

For the regime II bursts, we can set up an upper limit on t B using the light curves. 
Typically this limit is very low (say, < 10~ 4 ) since the afterglow light curves usually extend 
to very late epochs (say, > 1 day). The efficiency estimated for these bursts have even 
larger uncertainties because Ek depends on many unknown parameters, e.g. e e , e#, n, etc. 
Nonetheless, we perform some rough estimates with the following procedure. 

To keep v c above u x for a long period of time (say, days), one must have a small e# 
and/or a small n. Sometimes one even needs to lower e e if one is reluctant to lower n. A 
lower tB and a lower e e would lead to a larger Ek and hence, a lower r] 1 . Since we do not 
know which parameter is in operation, we fix e e = 0.1 and n = 0.1, and search e# downwards 
to find the highest e# that allows v m < v x < v c to be satisfied. In Table 4 we present our 
search results for the regime II bursts. The efficiency value listed does not represent the 
real value. There is no obvious reasons to argue why e e should be 0.1, why n should be 
0.1, or why e# is not even lower. Nonetheless, the value roughly indicates how Ek and 
look when the regime II spectral condition is satisfied. The results suggest that e B is rather 
low, in most cases lower than 10~ 4 . The estimated E K is typically very high, and hence, 
the radiative efficiency is very low. This suggests that either these GRBs are intrinsically 
inefficient or the ambient densities of at least these GRBs are quite low, say n < 0.1 cm -3 . 

5. Conclusions and discussion 

We present a detailed analysis on the prompt gamma-ray emission and the early X-ray 
afterglow emission for a sample of 31 GRBs detected by Swift and derived their radiative 
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efficiencies. The sample includes both long and short GRBs and both normal GRBs and 
softer XRFs. This allows us to investigate how the GRB radiative efficiency vary globally 
within different populations. We summarize our findings in the following. 

f. Due to the intrinsic limitation of the BAT instrument, it is very difficult to derive 
spectral parameters of the prompt emission and to extrapolate the BAT fluence to a broader 
band. We have developed a method by making use the hardness ratio information to derive 
spectral parameters assuming a Band function spectrum for all the bursts. This is probably 
the best effort with the available BAT data. The involved uncertainties is still large. It is 
also difficult to set errors to the quantities. The errors of r\ and of the Band function 
are not accounted for in our analysis, but those of E are included whenever possible. These 
errors are included to derive the errors of S y , which also include the uncertainties in the 
observations. Due to the intrinsic limitation of the method, the errors of S y are likely under- 
estimated. 

2. We compare our prompt emission data with those of BATSE (Preece et al. 2000) and 
HETE-2 (Lamb et al. 2005; Sakamoto et al. 2005) and generally find consistent results. Swift 
extends the fluence to the fainter regime and the hardness ratio to the softer regime with 
respect to BATSE. The derived properties of XRFs (including the 2:1 relative population 
between GRBs and XRFs) are generally consistent with those derived from the HETE-2 
data. 

3. The shallow decay component commonly detected in X-ray afterglows complicates 
the efficiency study. Previous analyses using late-time X-ray afterglow data (say, 10 hr) 
inevitably over-estimated Ek and hence, under-estimated the GRB efficiency if the shallow 
decay is due to energy injection. We define two characteristic time epoches, the putative 
fireball deceleration time (tdec) and the epoch when the shallow decay is over to study 
the efficiency problem. The efficiency derived at the former epoch [i.e. ^ 7 (tdec)] is likely 
the true efficiency for the models that interpret the shallow decay phase as due to energy 
injection (Zhang et al. 2006; Nousek et al. 2006; Panaitescu et al. 2006), while that derived 
at the later epoch [i.e. ?7 7 (i&)] is the true efficiency for some models that invoke precursor 
injection (Ioka et al. 2006) or delayed energy transition from the fireball to the circumburst 
medium (Kobayashi & Zhang 2006; Zhang & Kobayashi 2005). 

4. We investigate both the observationally defined gamma-to-X fluence ratio (Rj/x) 
and the theoretically defined efficiency (rjj) at the two epochs. The former invokes direct 
observables and therefore is subject to less uncertainties. The later involve theoretical mod- 
eling, which depends on the uncertainties of many unknown microphysics parameters. The 
error of the latter is therefore difficult to quantify. Our approach is to fix model parame- 
ters, and only derive the errors introduced from the uncertainties of observations and data 
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analyses. The errors of calculated efficiencies are therefore under-estimated. Nonetheless, 
some interesting features emerge from both the R-y/x analysis and the r/ 7 analysis. At t^, a 
shallow correlation between Ry/x or V-y an d the hardness ratio is found. This is consistent 
with previous findings (Soderberg et al. 2004; Lloyd-Ronning & Zhang 2004; Lamb et al. 
2005) that XRFs appear to be less efficient. However, the shallow correlation disappears 
when the analysis is performed at td ec - This suggests that if the shallow decay is indeed due 
to energy injection, XRFs have similar radiative efficiencies as normal GRBs. The apparent 
low efficiency inferred using late time data must be attributed to some other reasons. 

The result has important implications for understanding the nature of XRFs. In par- 
ticular, it disfavors the model that interprets XRFs as events similar to GRBs but having 
smaller Lorentz factor contrasts and therefore lower radiative efficiencies (e.g. Barraud et 
al. 2005), if the shallow decay phase is due to energy injection. It suggests that XRFs are 
dim because the initial kinetic energy is also low, and that they gain larger kinetic energies 
later through energy injection. While there is no straightforward reason for energy injection 
in the radial direction (e.g. due to a long-lived central engine or pile-up of slow ejecta) for 
XRFs, some geometric models of XRFs indeed expect energy "injection" from the horizontal 
directions. These models, such as the quasi-universal Gaussian-like structured jet model 
(Zhang et al. 2004a; Dai & Zhang 2005; Yamazaki et al. 2004), invoke relativistic ejecta 
with variable luminosities and possibly variable Lorentz factors in a wide range of angles. 
While GRBs correspond to the cases when observers view the bright core component (e.g. 
the line of sight is inside the typical Gaussian angle), XRFs are those cases when observers 
view the off-axis jet at a larger viewing angle. As long as the bulk Lorentz factor of the 
ejecta is still relativistic in these directions, the prompt emission of XRFs is weaker (because 
of the low energy in the direction) but with a comparable efficiency as in the core direction, 
since the initial kinetic energy in the direction is also low. Later as the jet is decelerated, the 
observer in the off-axis direction would progressively receive the energy contribution from 
the energetic core so that the light curve shows an early shallow decay (e.g. Kumar & Granot 
2003; Salmonson 2003). At later times, the effective kinetic energy in the viewing direction 
is enhanced, leading to an apparently high kinetic energy (and hence, an apparently low ra- 
diative efficiency). Such a picture might be consistent with the data. In fact, the long-term 
X-ray afterglow lightcurve of XRF 050416A could be modeled by such a model (Mangano 
et al. 2006). 

Other geometric models of XRFs have been also discussed in the literature. The model 
involving sharp-edge jets viewed off beam (Yamazaki et al. 2002) predicts an initial rising 
light curve and a steep decay after reaching the peak due to sideways expansion, which is 
not favored by the data. Granot et al. (2005) suggested a smoother edge, which in effect 
is similar to the Gaussian jet model (Zhang et al. 2004a). The two-component jet model 
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(Zhang et al. 2004b; Huang et al. 2004; Liang & Dai 2004; Peng et al. 2005) could be 
consistent with the data as long as the second component is not distinct enough to result 
in noticeable light curve features that are not detected. Power law structured jets (Zhang 
k Meszaros 2002b; Rossi et al. 2002) may also interpret XRFs (Jin & Wei 2004; D'Alessio 
et al. 2006), but the model predicts too large a XRF population (Lamb et al. 2005; Zhang 
et al. 2004a). Finally, the varying opening angle model for XRFs invoke very wide jets 
for XRFs (Lamb et al. 2005). However, there is no straightforward reason to expect the 
significant change of E K (and hence, 7/ 7 ) at early times in this model. 

It is possible that the shallow decay phase is not due to energy injection. Numerical 
simulations (Kobayashi & Zhang 2006) suggest that the time scale for a fireball to transfer 
energy from the ejecta to the medium is long. It could be that the observed shallow decay 
phase reflects the epochs during which the energy transfer is going on. A Poynting-flux- 
dominated flow may further extend the energy transfer period further due to the inability of 
tapping kinetic energy of the shell with the presence of a reverse shock (Zhang & Kobayashi 
2005). If this is the case, Exitb) is the relevant kinetic energy to define GRB efficiency, 
and XRFs are then indeed intrinsically inefficient GRBs (e.g. Barraud et al. 2005). More 
detailed early data and modeling are needed to reveal whether such a possibility holds. 

5. The absolute value of radiative efficiency is subject to the uncertainties of afterglow 
parameters. By inspecting the spectral index and the temporal decay index of the X-ray 
afterglows, we identify 22 bursts whose afterglow cooling frequency is below the X-ray band. 
In this regime the radiative efficiency is insensitive to parameters except for e e . Assuming 
e e = 0.1, we find that of most bursts are smaller than 10%, while ry 7 (tdec) ranges from 

a few percents to > 90%. Some bursts have a low efficiency throughout, and they correspond 
to those bursts whose X-ray afterglow light curve smoothly join the prompt emission light 
curve without a distinct steep decay component or an extended shallow decay component. 

The standard internal shock models predict that the GRB efficiency is only around 1% 
(Kumar 1999; Panaitescu et al. 1999). Our results indicate that some GRBs satisfy such 
a constraint. On the other hand, a group of GRBs with an early shallow decay component 
challenge the internal shock model if E K (t dcc ) is the relevant quantity to define the efficiency, 
as is required in most energy injection models. Suggestions have been made in the literature 
to increase the efficiency (e.g. Beloborodov 2000; Kobayashi & Sari 2001; Fan et al. 2004; 
Rees & Meszaros 2005; Pe'er et al. 2005). For those models that invoke Exitb) to define 
the efficiency (e.g. Ioka et al. 2006; Kobayashi & Zhang 2006), the data are consistent with 
the expectation of the internal shock model. The typical efficiency in such a case is 10% or 
lower. 



6. The three short GRBs (050724, 050813 and 051221A) have similar efficiencies as long 
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GRBs at both t dec and tb, despite their distinct progenitor systems (Gehrels et al. 2005; Fox 
et al. 2005; Villasenor et al. 2005; Hjorth et al. 2005; Barthelmy et al. 2005; Berger et 
al. 2005c; Bloom et al. 2006). Such a trend is being verified by more and more short GRB 
data, and will be reinforced by a larger sample of short GRBs in the future. 

7. Although most of the X-ray afterglows in our sample are above the cooling frequency, 
9 GRBs do have a cooling frequency higher than the X-ray band for a very long time, 
suggesting a very small e# and/or a low medium density n <C 0.1 cm~ 3 . An extremely low 
medium density may not be at odd for long GRBs. If some of these long GRBs explode in 
superbubbles created by precedented supernovae or GRBs, the ambient medium could reach 
such a low density. 

We thank an anonymous referee for helpful comments. This work is supported by NASA 
through grants NNG05GC22G, and NNG06GH62G (for BZ, EL) and by the National Natural 
Science Foundation of China (for EL, No. 10463001). 
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Table 1. Prompt Gamma- Ray data for the GRBs in our sample 
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Rcdshift references: (zl)Berger et al. (2005a), (z2) Bcrgcr ct al. (2005b), (z3)Fynbo ct al. (2005a), (z4)Fynbo et al. (2005b), 
(z5)Ccnko ct al. (2005), (z6)Foley et al. (2005), (z7)Berger et al. (2005c), (z8)Fynbo ct al. (2005c), (z9)Lcdoux ct al. (2005), 

(zlO) Berger & Soderberg (2005). 



a GRB duration in 15-150 keV. 

b The power law index and the reduced x 2 of the best fit to the BAT data. 

c The Hardness ratio is calculated by gamma-ray fluence in 50-100 keV band to that in 25-50 keV band. 
d The logarithm of the observed Gamma-ray fluence and its error in the 15-150 keV band. 

c The spectral parameters derived from the best Band- function fit with the constraint of HR° bs = HR mod , except for those bursts 
with marks. The errors of E v (in 90% confidence level) are derived from the best fits with Xspec package by fixing both Ti and Y2 ■ 

{ Logarithm of extrapolated gamma-ray fluence in 1 — 10 4 keV band. 

g The spectral parameters are taken from the Konus- Wind data (Golenetskii et al 2005a-d; Krimm ct al. 2006) 

h Most of the prompt emission was in one peak with 0.25s and therefore the burst was considered as a short burst(Barthelmy et al. 
2005). 



-41 - 



Table 2. Observations and Fitting Results of X-ray Afterglows 
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-6.64( 


.18) 


-7.24( 


.12) 


-6.76( 


.12) 


GRB050401 


.98( .05) 


.68( .05) 


1.35( .05) 


3.74( 


.04) 


-7.01( .10) 


-6.42( 


.11) 


-6.44( 


.11) 


-6.65( 


.12) 


GRB050406 


1.37( .25) 


.51( .10) 


- 


- 




-10.11( .33) 


-8.05( 


.11) 


-9.20( 


.14) 


-8.71( 


.06) 


GRB050416A 


.80( .29) 


.40( .05) 


.95( .05) 


3.18( 


.03) 


-8.51( .10) 


-7.68( 


.10) 


-7.61( 


.10) 


-7.55( 


.12) 


GRB050422 


2.23(1.09) 


.86( .04) 


- 


- 




-8.60( .20) 


-8.20( 


.16) 


-8.35( 


.17) 


-8.21( 


.16) 


GRB050502B 


.81( .28) 


.86( .01) 


- 






-8.38( .04) 


-7.97( 


.03) 


-8.12( 


.03) 


-7.98( 


.03) 


GRB050525 


1.07( .02) 


1.10( .05) 


1.45( .05) 


3.44( 


.05) 


-6.29( .11) 


-6.52( 


.12) 


-6.56( 


.12) 


-6.96( 


.13) 


GRB050607 


.77( .48) 


1.01( .05) 


2.46(1.65) 


5.11( 


.24) 


-7.97( .36) 


-8.06( 


.49) 


-7.99( 


.36) 


-8.00( 


.36) 


GRB050712 


.90( .06) 


.72( .01) 








-7.86( .17) 


-7.48( 


.14) 


-7.62( 


.15) 


-7.49( 


.14) 


GRB050713A 


1.30( .07) 


.65( .10) 


1.15( .10) 


4.00( 


.04) 


-7.55( .27) 


-6.86( 


.27) 


-6.96( 


.27) 


-6.89( 


.28) 


GRB050713B 


.70( .11) 


.30( .11) 


.97( .07) 


4.21( 


.23) 


-8.23( .20) 


-6.66( 


.24) 


-7.06( 


.20) 


-6.59( 


.29) 


GRB050714B 


4.50( .70) 


.50( .06) 








-8.17( .32) 


-6.82( 


.26) 


-7.36( 


.27) 


-6.86( 


.26) 


GRB050716 


.33( .03) 


1.06( .07) 








-7.41( .38) 


-7.58( 


.31) 


-7.51( 


.32) 


-7.57( 


.31) 


GRB050717 


1.15( .10) 


1.49( .01) 








-6.63( .05) 


-8.00( 


.04) 


-7.47( 


.04) 


-7.96( 


.04) 


GRB050721 


.74( .15) 


1.28( .01) 








-6.78( .06) 


-7.60( 


.05) 


-7.30( 


.05) 


-7.58( 


.05) 


GRB050724 


.95( .07) 


.90( .0.1) 








-9.21( .33) 


-7.86( 


.07) 


-8.56( 


.14) 


-8.21( 


.06) 


GRB050726 


.94( .07) 


.91( .10) 


2.45( .10) 


4.06( 


.04) 


-7.45( .11) 


-7.30( 


.12) 


-7.28( 


.11) 


-7.96( 


.14) 


GRB050801 


.72( .54) 


1.10( .02) 








-7.70( .09) 


-7.99( 


.06) 


-7.89( 


.07) 


-7.99( 


.06) 


GRB050802 


.91( .19) 


.65( .10) 


1.65( .10) 


3.76( 


.04) 


-7.46( .10) 


-6.80( 


.11) 


-6.82( 


.10) 


-7.26( 


.14) 


GRB050803 


.71( .16) 


.61( .10) 








-7.58( .28) 


-6.77( 


.08) 


-7.04( 


.14) 


-6.65( 


.06) 


GRB050813 


2.42( .89) 


.70( .10) 


2.12( .10) 


2.26( 


.12) 


-9.40( .13) 


-9.29( 


.20) 


-10.69( .30) 


-11.81( .36) 


GRB050814 


1.08( .08) 


.62( .06) 








-8.63( .34) 


-7.51( 


.28) 


-7.97( 


.29) 


-7.59( 


.28) 


GRB050819 


1.18( .23) 


-.06( .06) 


.65( .10) 


4.26( 


.16) 


-10.85( .10) 


-8.20( 


.11) 


-8.89( 


.10) 


-8.04( 


.14) 


GRB050820A 


.87( .09) 


1.18( .01) 








-6.07( .03) 


-6.46( 


.01) 


-6.27( 


.02) 


-6.45( 


.01) 


GRB050822 


1.60( .06) 


.46( .11) 


.95( .04) 


4.14( 


.31) 


-8.20( .29) 


-7.11( 


.33) 


-7.37( 


.29) 


-7.03( 


.38) 


GRB050826 


1.27( .47) 


1.10( .01) 








-7.41( .05) 


-7.70( 


.04) 


-7.60( 


.05) 


-7.70( 


.04) 


GRB051221A 


1.04( .20) 


1.07( .03) 








-7.99( .11) 


-8.13( 


.06) 


-8.12( 


.06) 


-8.19( 


.05) 



a X-ray spectral index 

h a\ and ai are the temporal decay indices before and after the break time (t;,). If a light curve is fitted by a simple power law only 
a\ is available 

c The values of S x at a given time are calculated by the flux times the corresponding time, in units of erg cm~ 2 . 
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Table 3. The derived kinetic energies and gamma-ray Efficiencies at td ec and tb for the GRBs in 
the spectral regime v x > max(z/ m , u c ) by assuming (e e , e^) = (0.1, O.Ol)^ and 

(e e ,e B ) = (0.1,0.0001)( 2 ) 



GRB log (tdec) ^(tdec)' logB^(t6) ^fo)* log^ 2) (t dcc ) ^ 2) (t dcc )* logfi^fo) ^ 2) (t 6 )* 







(erg) 


(100%) 




(erg) 


(100%) 




(erg) 


(100%) 




(erg) 


(100%) 


GRB050126 


51 


86( 


.84) 


57.78(48.81) 


53 


60 ( 


.48) 


2.39( 2.84) 


52 


26 ( 


.84) 


35.27(45.68) 


54 


00 ( 


.48) 


.97(1.16) 


GRB050128 


54 


48 ( 


.16) 


4.03(1.51) 


55 


01( 


.17) 


1.22(.50) 


54 


49 ( 


.16) 


3.99(1.49) 


55 


02( 


.17) 


1.21(.50) 


GRB050219A 


52 


84( 


.39) 


57.21(22.20) 


54 


05 ( 


.30) 


7.57( 4.85) 


52 


85 ( 


.39) 


56.09(22.33) 


54 


07( 


.30) 


7.25(4.66) 


GRB050315 


51 


69( 


.06) 


93.55(.96) 


54 


59( 


.08) 


1.79(.33) 


52 


09 ( 


.06) 


85.24( 2.00) 


54 


99( 


.08) 


.72(.13) 


GRB050319 


52 


45( 


.10) 


~ 70 


55 


59( 


.14) 


~ 0.2 


52 


85 ( 


.10) 


~48 


55 


99( 


.14) 


~.l 


GRB050401 


55 


05 ( 


.10) 


2.79(.70) 


55 


66 ( 


.11) 


.70(.19) 


55 


06 ( 


.10) 


2.76(.69) 


55 


66 ( 


.11) 


.70(.18) 


GRB050406 


50 


67( 


.28) 


~80 


54 


38 ( 


.09) 


~0.1 


50 


99 ( 


.28) 


~66 


54 


69 ( 


.09) 


~.04 


GRB050416A 


52 


34( 


.10) 


~5 


53 


18( 


.10) 


~.8 


52 


35 ( 


.10) 


~5 


53 


18( 


.10) 


~0.8 


GRB050422 


51 


80( 


.16) 


68.99(12.56) 


53 


83 ( 


.13) 


2.01( 1.06) 


52 


20( 


.16) 


46.97(14.62) 


54 


23( 


.13) 


.81(.43) 


GRB050502B 


53 


40 ( 


.04) 


3.99(.99) 


53 


84( 


.03) 


1.51(.37) 


53 


41 ( 


.04) 


3.95(.98) 


53 


84( 


.03) 


1.50(.37) 


GRB050525 


53 


50 ( 


.11) 


7.64(1.78) 


53 


45 ( 


.11) 


8.49(2.03) 


53 


56 ( 


.11) 


6.61( 1.55) 


53 


52( 


.11) 


7.36( 1.78) 


GRB050607 


53 


82( 


.36) 


4.11(4.06) 


53 


75 ( 


.49) 


4.78(5.86) 


53 


82( 


.36) 


4.06( 4.02) 


53 


75 ( 


.49) 


4.72( 5.80) 


GRB050712 


54 


18( 


.14) 


1.63(.86) 


55 


03 ( 


.14) 


.24(.13) 


54 


19 ( 


.14) 


1.61(.85) 


55 


03 ( 


.14) 


.23(.12) 


GRB050713A 


52 


92( 


.24) 


68.02(12.41) 


54 


37( 


.24) 


7.07(3.76) 


53 


19( 


.24) 


53.84(14.18) 


54 


63 ( 


.24) 


4.00( 2.20) 


GRB050714B 


50 


80( 


.25) 


~96 


53 


49( 


.21) 


4 


51 


20( 


.25) 


~ 89 


53 


89 ( 


.21) 


~ 2 


GRB050717 


53 


91( 


.05) 


14.11( 6.38) 


53 


221 


.04) 


44.47(12.90) 


54 


05 ( 


.05) 


10.65(5.01) 


53 


36 ( 


.04) 


36.74(12.14) 


GRB050724 


50 


40 ( 


.29) 


69.23(14.68) 


50 


81( 


.07) 


46.54(5.16) 


50 


40 ( 


.29) 


68.98(14.75) 


50 


81( 


.07) 


46.26(5.15) 


GRB050813 


51 


03 ( 


.10) 


44.81(16.69) 


51 


45 ( 


.16) 


23.63(13.13) 


51 


43 ( 


.10) 


24.43(12.46) 


51 


85 ( 


.16) 


10.96(7.10) 


GRB050814 


52 


17( 


.32) 


67.78(16.47) 


53 


59( 


.27) 


7.38(4.31) 


52 


25 ( 


.32) 


63.80(17.42) 


53 


67( 


.27) 


6.26(3.70) 


GRB050819 


49 


95 ( 


.09) 


~ 99 


53 


02( 


.10) 


~7 


50 


H( 


.09) 


~98. 


53 


18( 


.10) 


~ 5 


GRB050822 


52 


43 ( 


.23) 


~69 


54 


59( 


.27) 


~ 1.55(.94) 


52 


83 ( 


.23) 


~47 


54 


99( 


.27) 


~ .6 


GRB050826 


53 


01( 


.05) 


> 12 


53 


80( 


.04) 


> 2 


53 


25 ( 


.05) 


> 7 


54 


04( 


.04) 


> 1 


GRB051221A 


51 


92( 


.11) 


32.56( 7.83) 


51 


90( 


.06) 


33.57(6.54) 


51 


96 ( 


.11) 


30.61( 7.57) 


51 


94( 


.06) 


31.59( 6.33) 



* The derived r\ is insensitive to the rcdshift. For those bursts whose rcdshifts are not available we take z = 2 in the 
calculation. The errors of r\ are calculated by considering only the errors from the extrapolated gamma-ray fluences and the 
observed X-ray fluxes. Therefore, the errors of r\ only reflect the observational errors. Ex is sensitive to microphysics 
parameters that are poorly constrained. The true errors of rj should be significantly larger than what are reported here. 
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Table 4. Derived kinetic energies and gamma-ray efficiencies 
at tdec and % for the GRBs in the spectral regime 
v m < v x < v c with e e = 0.1 



GRB 










Ek,52 


(tdec) 




Ek,52 


(*6) 














(er 


g) 


(100%) 


(erf 


;) 


(100%) 


GRB050215B 


3 


X 


io- 


-6 


53.65( 


1.49) 


~ 1 


54.17( 


1.40) 


~0.4 


GRB050713B 


4 


X 


10" 


-6 


52.01( 


.35) 


90.68( 7.19) 


54.98( 


.42) 


1.03( 1.02) 


GRB050716 


3 


X 


icr 


-7 


55.96( 


.76) 


.20( .36) 


55.09( 


.62) 


1.48( 2.10) 


GRB050721 


9 


X 


io- 


-5 


54.55( 


.10) 


2.11( .66) 


53.54( 


.08) 


18.29( 4.35) 


GRB050726 


5 


X 


io- 


-5 


53.22( 


.16) 


>48 


54.38( 


.18) 


>6 


GRB050801 


2 


X 


io- 


-4 


52.84( 


.15) 


~ 11 


52.93( 


.10) 


~9 


GRB050802 


2 


X 


io- 


-5 


52.95( 


.16) 


29.60( 9.72) 


54.79( 


.17) 


.61( .29) 


GRB050803 


6 


X 


io- 


-6 


53.12( 


.48) 


29.37( 23.78) 


54.79( 


.14) 


.88( .37) 


GRB050820A 


7 


X 


io- 


-7 


56.50( 


.04) 


.14( .04) 


56.22( 


.02) 


,27( .07) 



